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Abstract. The motion of a fully ionized plasma of electrons and ions is gen- 
erally governed by the Vlasov-Maxwcll-Landau system. We prove the global 
existence of solutions near Maxwellians to the Cauchy problem of the sys- 
tem for the long-range collision kernel of soft potentials, particularly including 
the classical Coulomb collision, provided that initial data is smooth enough 
and decays in velocity variable fast enough. As a byproduct, the convergence 
rates of solutions are also obtained. The proof is based on the energy method 
through designing a new temporal energy norm to capture different features 
of this complex system such as dispersion of the macro component in R 3 , sin- 
gularity of the long-range collisions and regularity-loss of the electromagnetic 
field. 



1. Introduction 

1.1. Kinetic equations for a fully ionized plasma. The motion of a fully 
ionized plasma consisting of only two species particles (electrons and ions) under 
the influence of the self-consistent Lorentz force and binary collisions is governed 
by the kinetic transport equations 

d t F+ + £ • V X F+ + (E + dxB)- V e F+ = Q(F+, F+) + Q(F+, F_), 

dtF- + £ • V K F_ - (E + £ x B) ■ V C F_ = Q(F_ , F+) + Q(F_, F_), 

coupled with the Maxwell equations 

d t E - V x x B = - [ i(F + - F_) d£, 

d t B + V, x E = 0, (1.2) 
V x -E= ( (F+-F-)d£, V x -B = 0. 



(1.1) 



Here F± = F±(t, x, £) > stands for the number densities of ions (+) and electrons 
(— ) which have position x = (xi, £2,2:3) € K 3 and velocity £ = (£17^2,^3) € ^ 3 a t 
time t > 0, and E(t, x), B(t, x) denote the electro and magnetic fields, respectively. 
The initial data of the system of equations (|l.ll) - (|1.2[) is given by 

F ± (0,x,0=F 0t± (x,0, E(0,x)=E o (x), B(0,x) = B (x), 

satisfying the compatibility conditions 

V x -E = ( (F ,+ - F ,_) d£, V, • B Q = 0. 

it 3 

Notice that for simplicity of presentation and without loss of generality, all the 
physical parameters appearing in the system, such as the particle masses and the 
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light speed, and other involving constants, have been chosen to be unit. Moreover, 
the relativistic effects are neglectd in our discussion of the model. This is normally 
justified if the plasma temperature is much lower than the electron rest mass, cf. |13[ 
Section 3.6]. We refer interested readers to p~3l Chapters 2 and 3] or [16j Chapter 6] 
for the non dimensional representation of the system (jl.ip - (|1.2[) and also its physical 
importance in the study of plasma transport phenomena. 

The quadratically nonlinear operator Q in (jl.lj) is described by the Landau 
collision mechanism between particles in the form of 

Q{F,G)=V V { f <!>(£- t*)[VeF(OG&)-F(t)Vs,G(Z*)]d!;* 
Ur 3 

The non-negative matrix <p denotes the Landau collision kernel 

4> ij (0 = CM 1+2 {% - ||} , -3 < 7 < -2, 

for a constant > 0. Note that 7 = —3 corresponds to the Coulomb potential 
for the classical Landau operator which is originally found by Landau (1936). For 
the mathematical discussions of the general collision kernel with 7 > —3, see the 
review paper [26) by Villani. 

From now on, for simplicity of presentation, it is convenient to call (|1.1[) - (|1.2I) 
the Vlasov-Maxwell-Landau system. In this paper we aim at proving the global 
existence of solutions to the Cauchy problem of the system near a global Maxwellian 
under some conditions on initial data. 

1.2. Reformulation. Write the normalized global Maxwellian as 

M = Me) = (2.)- 3 / 2 e-l«l 2 / 2j 

and set the perturbation in the standard way 

Use [•,•] to denote the column vector. Set F = [F+,F—] and / = [/+,/_]. Then 
the Cauchy problem (|l.lj) - ([1.2[) can be reformulated as 

= |£.£/ + r(/,/), 

<d t E-V x xB = -( £m 1/2 (/+-/-K, (1.3) 
Jm. 3 

d t B + V x x E = 0, 

V x -E= f M 1 / 2 (/ + -/_)dC, V, -B = 0, 
y jr3 

with initial data 

/±(0,z,0 = /o,±M, E(0,x)=E (x), B(0,x) = B o (x), (1.4) 
satisfying the compatibility condition 

V x • E a = f ^ 1/2 (/o,+ - fo,-) df, V, ■ B = 0. (1.5) 

Js. 3 
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Here, qo = diag(l,— 1), qi = [1, — 1], and the linearized collision term Lf and the 
nonlinear collision term r(/, /) are respectively defined by 

Lf=[L+f,L-f\, T(f,g) = \T + (f,g),T-{f,g)], 

with 

L±f = -2 M - 1/2 0(M 1/2 /±, /x)-/i- 1/2 Q(/i, M V2 {/± + /*}), 

r ± (/,.9) = M- 1/2 Q(M 1/2 / ± ,M 1/2 3 ± ) + M- 1/2 Q(M 1/2 /±,M 1/2 3 T ). 

1.3. Macro projection, weights and norms. As in [TTJH], the null space of 
the linearized operator L is given by 

N = span {[1,01m 1 / 2 , [0, l}^ 1/2 , [tiM^ 2 (1 < « < 3), [|£| 2 , |£| V /2 } ■ 

Let P be the orthogonal projection from L 2 x L 2 to N. Given /(i, ai,£), one can 
write P as 

P/ = a + (t,x)[l,0]fj} /2 +a-(t,x)[0, l]/i 1/2 

+ ^ b t (t,x)[l, l]^ 2 + c(t, x)[l, - 3)^ 2 , 

i=l 

where the coefficient functions are determined by / in the way that 
a± = (n 1/2 J±) = (v 1/2 ,P±f), 

h = \{tm 1/2 J+ + f-) = (tm 1/2 ,P±f), 

c = 3) M V2 ) /+ + ; _) = l m * 3) M 1/2 ,P ± />. 

In what follows, we introduce the weight functions and norms used for the pre- 
sentation of the main result later on. Fix a constant i? with < i? < 1/4. Define 

w r ,x = w T , x (t,0 = (0 (7+2)T ex P { (Y^(0 2 } , 

where constants r £ R and A > are two parameters which may vary in different 
places. For / = f(t,x,£), define 



\f(x)\ 2 T ^= w 2 TtX (t,0\f\ 2 dt, ||/|| 2 , A =/ \f(x)\ 2 !X dx, 

JR 3 JR 3 

and 

\m\h, T ,x = E J R3 <x(t,o {^w/+^||i/i 2 } 

ll/llD.r.A= / l/WllJ.r.A^, 

JR 3 

where = % denotes the velocity derivative with respect to and cr^ = cr lJ '(£) 
is the Landau collision frequency given by 

cr« (0 = pi * = [ <t> ij (£ -e«>K.) d £* ■ 

JR 3 
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To study the global existence through the energy method, the temporal energy 
functional and the corresponding dissipation rate are defined by 

6v A a(*)~ E \\tfW)\\M-i,x + ME,B)\\* BJt , (1.6) 

\a\ + \0\<N 



and 



T> N ,t,x(t)= E ll^U-n/Wllo, 

|a|+l/3|<JV 

+ E l|V.^(a ± ,6, C )|| 2 



|a|<iV-l 

+||a+ - a_f + H^ll^-i + ||V*fl||^_, 

+ fl E iiwu-n/wiif^A, (i-7) 

where the integer iV > and £ > are parameters which may differ in different 
places. For simplicity, we write w T = uv.o for A = 0, and £jv(t) = £jv,o,o(i) for 
^ = A = 0, and likewise for others. Moreover, regarding £N,i,\{t) and T>N,£,\{t), 
whenever the case A > occurs, we always require £ — N > 0; it means that the 
algebraic weight factor is always of positive power when the exponential weight 
factor is present. Notice that the last term of T> N ^\(t) in (|1.7I) disappears when 
A = 0. 

Let constants Nq and £q be fixed properly large, and let constants Ao > and 
eo > be fixed properly small; the choice of Nq, £q, \q and eo can be seen in the 
late proof. Set N\ = |iVo and £\ = \t§. The temporal energy norm X{t) is defined 

by 

X{t) = sup {£ Nl (s) + (1 + s)h Nl - 2 (s)} 

0<s<t 

l + ep 

+ sup{(l + s) 2 £ Nu e u \ (s) + £jVi-ia,a (s) 

0<s<t 

+ (1 + s)l^jv 1 -3A-l.Ao( s )} 

3 

+ sup {£jVoA,A ( s ) + (1 + s )'^Afo,«o-l,A (s)} 
0<s<t 

+ Bup{(l + *) a ( 1 +*)||V a (S,B)( S )||^ _ 1 }. (1.8) 

0<s<t 

Notice that in order for the algebraic weight factor to gain the large enough positive 
power when there is the exponential weight factor exp{Ao(£) 2 /(l wc a l so let 

Co — 3Nq be properly large. 



1.4. Main result. The main result of the paper is stated as follows. 

Theorem 1.1. Assume —3 < 7 < —2. Take < < 1/4, and also take con- 
stants Nq, £q properly large with £q — 3Nq properly large, and constants Ao > 0, 
eo > properly small. Fix a constant £2 > § + Let fo = [/o,+ ,/o,~] satisfy 
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F±(0, x, = + /i 1 ^ (0/o,± (.t, > 0. If 

Yo= ^ ll^/olll/il-^Ao + ^2 ll^/o|l|/3|-fi,A 
\a\ + \f3\<N |a| + |/3|<iVi 

+ \\(E ,B )\\ H N inZl + \\W-eJo\\ Zl (1-9) 

is sufficiently small, then there are appropriately defined energy functionals £N.£.\{t) 
appearing in X(t) such that the Cauchy problem (|1.3p . (|1.4|) . (| 1 . 5j) of the Vlasov- 
Maxw ell-Landau system admits a unique global solution (f{t,x,^),E{t,x),B{t,xy\ 
satisfying F±{t,x,£) = + fJ, 1/2 (Of±{t, x,£) > and 

X{t) < Y 2 , (1.10) 

for all time t > 0., 

Remark 1.1. The following are several points to remark on this theorem: 

(a) Although only the soft potential case —3 < 7 < — 2 is considered here, the 
hard case for 7 > — 2 could be much simpler. 

(b) Even when 1? is strictly less than 1/4, it can sitll be shown from the late 
proof that 

£ (ir/u + \\d«(E,B)\\) 

l<\a\<N a 

decays in time with the rate (1 + i)~ 5 / 4 , cf. (|4.34[) . 

(c) By the definition of X(t), the uniform-in-time inequality (|1.10[) implies that 
the weighted high-order energy functional £jvi,£i,a (0> particularly 

E w-nfiif,,-^, 

|a| + |j8|=JVi 

may increase in time with the rate (1 + i)( 1 + e °)/ 2 . 

(d) S'ettmg iVi = |iV anrf £1 = *s jwsf /or simplicity of presentation. The 
general choice of N\, t\ in terms of No, £q is possible. In addition, for 
brevity, Nq, £q and £q — 3iVo are assumed to be properly large. We would 
not track in this paper the critical values of Nq, £q and all other param- 
eters. However, it should remain an interesting problem to design a new 
energy norm with the optimal choice of regularity and velocity integrability 
on initial data in order to ensure the global existence of the Cauchy problem. 

(c) The similar approach developed in this paper could be immediately applied 
to the case of the periodic box. In that case, we need to assume all the con- 
servation laws of the system so that the Poincare inequality can be applied 
to deal with the zero-order dissipation of the macro component Pf and the 
magnetic field B. 

In what follows we mention some work only related to this paper; interested 
readers may refer to them for more references therein. Notice that there are different 
approaches in establishing the mathematical theories on the Landau equation, see 
[D 01 [T7J [55]. In the perturbation framework, Guo [TU] firstly established 
the global existence for the purely Landau equation with Coulomb potentials in 
the absence of any force; see also [15]. Very recently the same author [11] made a 
further progress for the Vlasov-Poisson-Landau system on the periodic box when 
the self-consistent potential force is present; see [7] and [24] for generalizations of 
the result to the case of the whole space. We pointed out that Duan- Yang- Zhao [7] 
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used a different approach arising from the study of the Vlasov-Poisson-Boltzmann 
system 0|6]. 

When the plasma transport is effected by the Lorentz force coupled with the 
Maxwell equations, there are two cases in which the global existence of solutions 
near the global Maxwellian have been well studied. One case is to take the Landau 
operator with 7 > — 1. In this case the problem was solved by Guo [12] although 
the Vlasov-Maxwcll-Boltzmann system for only the hard-sphere model is considered 
there. Notice from Lemma [2.11 that if 7 > — 1, i.e. 7 + 2 > 1, then compared to 
the soft potential case, the linearized Landau operator has the stronger dissipative 
property which can control those nonlinear terms with the velocity-growth rate 
|£|, typically occurring to E ■ £/. The other case is the relativistic version of the 
Vlasov-Maxwell-Landau system studied by Strain-Guo [23]. In this case, due to 
the boundedness of the relativistic velocity and the special form of the relativistic 
Maxwellian, the velocity growth phenomenon in the nonlinear term disappears, 
and instead the more complex property of the relativistic collision operator was 
analyzed there. 

Though there are a few results mentioned above, it still remains unknown to 
obtain the global existence for the Vlasov-Maxwell-Landau system for the Coulomb 
potential in the classical sense because of the quite complex property of the system 
that we will point out in more detail next subsection. To the knowledge of our 
best, Theorem 1 1.1 1 is the first result in this direction. The proof of Theorem 11.11 is 
based on the energy method, cf. [TUl [12] or [H], together with a new quite delicate 
bootstrap argument. In particular, the introduction of the temporal energy norm 
X(t), which is the main strategy of the proof, captures most of important features 
of the coupling system. 

1.5. Difficulty and idea in the proof. The Vlasov-Maxwell-Landau system un- 
der consideration has the following three typical features which result to different 
mathematical difficulties: 

• the degeneration of dissipation at large velocity for the linearized Landau 
operator with soft potentials; 

• the velocity-growth of the nonlinear term, as mentioned before; 

• the regularity-loss of the electromagnetic field. 

The first feature makes it impossible to control the velocity derivative of the linear 
transport term £ • V x f without any velocity weight. To deal with it, the algebraic 
weight factor (£)(t+ 2 )I^I depending on the order of velocity differentiation was in- 
troduced in [10] . This, however, induces another difficulty when estimating the 
nonlinear term (E + £ x B) ■ Vg /, since the term contains one velocity differen- 
tiation so that the extra velocity-growth with rate |£|~(i'+ 2 ) is produced. Notice 
that this new trouble as well as the obvious velocity-growth in E ■ £/ happen to 
the nonlinear term only. [5] then introduced a time-velocity dependent exponential 
weight factor exp{Ao(£) 2 / (1 +t)®} which from the weighted estimate on dtf indeed 
leads to the following additional good term in the dissipation rate 



See also 0[7]. Thus, as long as those nonlinear velocity-growth terms contains a 
portion decaying in time faster than (1 +t) 1+ ®, they can be controlled by the above 
dissipation term. 



j J "<£>Hm,Ao(*.OIW - P}f\ 2 dxdZ 
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The third feature mentioned before can be seen from the definition of the dissipa- 
tion rate functional T>N,e,x{t) in which the L 2 norm of A^th-order spatial derivatives 
of (E, B) is missing. Here, the regularity-loss results essentially from the coupling 
of the hyperbolic Maxwell equations but not from the technique of the approach; 
see [4] for the analysis of the Green's function of the damping Eulcr-Maxwell sys- 
tem. Due to the regularity-loss, two more difficulties appear. One difficulty is that 
one can not expect derivatives of the electromagnetic field (E, B) of all orders up 
to the largest number N\ to decay time fast enough, so that the estimate on those 
nonlinear velocity-growth terms is still a problem. To solve it, when estimating the 
weighted inner product term 

E wl E -a-(E+z*B). v e /}, w 2 m _ iuXo d$f), 

M+|/8|<JVi 

we use the time-decay property for the only low-order derivatives of both (E, B) and 
/. Notice that compared to the high-order energy functional £Ni,ti,\o(t) 01 fi tnc 
low-order one must be assigned with the higher velocity weight to absorb the extra 
velocity-growth factor. That is the reason why we introduce into the X(t) norm 
two energy functionals £ jVi,ii,A (*) an d £N ,e ,\ (t) with the approximate choice of 
N and I. 

The second difficulty due to regularity-loss is the control of inner product terms 

|a|=JVi 

which arises from the weighted estimate on d a f. To deal with it, we use the time- 
weighted energy estimate with the time rate of negative power; the similar technique 
has been used in In fact, starting from the Lyapunov inequality 

j£ Nl (t) + KV Nl {t) < h.o.t, 
where h.o.t. denotes the high order terms only, it follows that 
d 



di [(l + t)- eo £ Nl (t)}+ K (l+t)- e "V Nl (t) 



(i + t) 1+eo 

Therefore, provided that the term on the right is time integrable, one can recover 
the following new dissipation term 

V \\d a (E,B)\\ 2 . 

f1 4- t) 1+e o ^ " v 

Although this dissipation term is degenerate in large time, it is enough to control 
other trouble terms related to the A^th-order derivatives of (E, B). 

We finally mention Theorem 13.11 concerning the time-decay property of the lin- 
earized system. Theorem l3.1l not only plays a key role of dealing with the dispersion 
of the macro component Pf in M 3 due to the degeneration of L, but also its proof, 
particularly the Lyapunov inequality (|3.10j) . fully reveals the optimal dissipative 
structure of the linearized system, which further motivates the design of the energy 
norm X(t). We remark that inspired by (|3.10p . it would be an interesting problem 
to consider the spectrum [23 [2] or Green's function [19] of such complex system. 
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The rest of the paper is arranged as follows. In Section [21 we list some known 
facts for the macro structure of the system and also the basic estimates on L and 
r. In Section [3J we obtain the time-decay property of the linearized homogeneous 
system. In Section [4l we present series of lemmas for the a priori estimates on the 
solution and finish the proof of Theorem 11.11 at the end. 

1.6. Notations. Throughout this paper, C denotes some generic positive (gener- 
ally large) constant and k denotes some generic positive (generally small) constant, 
where both C and k may take different values in different places. A < B means 
that there is a generic constant C > such that A ^ CB. A ~ B means A < B and 
B < A. We use L 2 to denote the usual Hilbert spaces L 2 = L 2 ^ or L 2 with the norm 
|| • || , and use (•, •) to denote the inner product over L 2 ^ or L 2 . For q > 1, the mixed 
velocity-space Lebesgue space Z q = L 2 (L1) = L 2 (]R|; L q (M^,)) is used. For multi- 

indices a = (a l5 a 2 ,« 3 ) and f3 = fafafo), d a p = = d^d^d^d^d^d^ . 

The length of a is \a\ = a\ + a 2 + ct 3 and similar for \/3\. 



2. Preliminary 



2.1. Macro structure. Consider the following linearized Vlasov-Maxwcll-Landau 
system with a non- homogeneous source 5* = [S+(t, x, £),<5L(t, x, £)]: 



f d t f ± + £ • V s /± T E ■ C/i 1/2 + L±f = S±, 
d t E-V x xB = - f ^ 1/2 (/+ - /-) df, 

JR 3 



d t B + V x x E = 0, 



(2.1) 



p 1/a (/+-/-K> V x -B = 0. 



Taking velocity integrations of the first equation of (|2.1[) with respect to the velocity 
moments 

M 1 / 2 , ^ 1/2 ,* = 1,2,3, I(| C | 2 -3)M 1/2 , 

6 

one has 

d t a± + V, • b + V, • (^ 1/2 , {/± - P±}/) - (/ /2 , S±), 
dt[bi + (^M 1/2 , U± - P±}f)] + d l (a ± + 2c) T E< 
+V a • (^M 1/2 , U± - P±}/> = (&M 1/2 , S±-L ± f), 

9 t c+ l -m 2 -zV 2 Ai±-p±}f) +l^ x -b 
+\^ x ■ m\ 2 - 3)^ 1/2 , {/± - p±}f) = line 3)^, s ± -L ± f), 

where d l = d Xi denotes the spatial derivative with respect to x,, and we have set 
/ = [J_|_, 7_] with I±f — f±- As in [51 H] , we define the high-order moment functions 
e(/±) = (e«(/ ± )) 3X 3 and A(/±) = (A 1 (/ ± ),A 2 (/±) ! A 3 (/ ± )) by 

e«(/±) = - i)m 1/2 ,/±>, Ai(/±) = - 5)6m 1/2 ,/ ± ). 
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Further taking velocity integrations of the first equation of ()2.1|) with respect to 
the above high-order moments one has 

d t [Q u ({I± - P±}f) + 2c] + 2d% = 0«(r± + S±), 
d t B l3 ({I± - P±}f) + d j h + d% + V, • (^ 1/2 , {I± - P±}f) 

d t Ai({I± - P±}f) + d l c = Ai(r± + S±), 

where r± = -£ • V K {I ± - P±}f - L±f. 

In particular, for the nonlinear system (|1.3p . the non- homogeneous source S = 
[S+(t,x,£), S-(t,x,£)] takes the form of 

S± = ±\e -U±T(E + t*B)- v c /± + r±(/, /). 
Then, it is straightforward to compute from integration by parts that 

(m 1/2 ,s±) =0, 

(^ 1/2 , S±) = ±Ea ± ± b x B±(£/i 1/2 , {J± - P±}/) x B 

+(^ 1/2 ,r±(/,/)), 

i((|e| 2 - 3)/i 1/a , 5±> = ±ift • E± {I± - P±}f) ■ E 

+(^(ICI 2 -3)/ /2 ,r ± (/,/)). 

2.2. Basic estimates on L and T. In this section, we state two lemmas about 
some basic properties of the Landau operator. Given a vector-valued function 
u = (ui, U2, U3), define 



lei J ki 



i ki M i^i ■ 



Concerning the equivalent characterization of the dissipation rate and the dissipa- 
tive property of the linearized Landau operator, one has the following lemma; see 
H Hi EH for the detailed proof. 

Lemma 2.1 ([ID]). If holds that 



I/Id.t.a 



(i + KI)*p*v c / 



(i + iei)^u-mv«/ 



"Y + 2 



Moreover, L is a nonnegative definite self-adjoint operator, and there exists k > 
smc/i that 

(LfJ)> K \{I-P}f\ 2 u . 
The following lemma states the weighted estimate on Lf and /). 
Lemma 2.2 ([22]). There are n > and C > swc/i t/iaf 

(X/, < A (t, 0/> > K\f\h, T ,X C |X{|f |<2C7}/|! ■ 
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Moreover, let |/3| > 0, r 
is C-q > such that 



£ with I > 0. Then, for 77 > small enough, there 



l^'l=l/3| 

-C?) E I^(9'/Id,I/3'|-^A- 

,4nd a/so, /or TV > 8, |a| + \/3\ < N and t = \P\-£ with £ > 0, it holds that 



E 



< 



\a'\ + \P'\<N 
0"<P'<P 



D.r, A 



D,T 



T.X 



mi-. 



pJ 3 \n,T,\ 



(2.2) 



Notice that since the coefficient X/(l + t)^ in front of (£) 2 in the exponential 
part of the weight function w T ,\{t,£) is bounded uniformly in t > 0, the proof of 
the above lemma follows directly from the same argument used in |22) . 



3. Linearized analysis 



Consider the Cauchy problem on the linearized Vlasov-Maxwell-Landau system 
with a source S = S(t,x,£) = [S+(t, x,£),S- (t, x, £)]: 



f 9 t f + £ ■ Vje/ — -B ■ £M 1/2 <Zi + L/ = 5, 

w-v I xB=-(( (1 i / 2 ,/ + -n 

V x -£ = ( M 1/2 , /+-/-), V x -fl = 0, 
, (/>-Bj-B)k=o = (fo,Eo,B ), 

where initial data [fa, Eg, Bq\ satisfies the compatibility condition 
V x -E a = [ ^ 1/2 (/o.+ -/o.-)^, V x -B = 0, 
and the source term S is assumed to satisfy 

H 1/2 {S+- S_)df = 0. 



(3.1) 



(3.2) 



/ 



To consider the solution to the Cauchy problem (|3.1[) . for simplicity, we denote 
U = [/, E, B], Uo = [fo, Eq, Bo] so that one can formally write 

U(t) = A(t)U + [ A(t-s)[S(s), 0,0] ds, 



where A(t) is the linear solution operator for the Cauchy problem on the linearized 
homogeneous system corresponding to p.l[) in the case when 5 = 0. 
For the linearized homogeneous system, we have the following result. 
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Theorem 3.1. Let S = 0, and let [f,E,B] be the solution to the Cauchy prob- 
lem (|3.1|) . (|3.2|) of the linearized homogeneous system. Define the velocity weight 
function w = by 

w(a) = (0- 2 ¥ (3-3) 
Then, for £ > and a > with m = \a\, 

\\w e d a f\\ + \\d a {E,B)\\ 

<(l + t)-^(\\w e+e '° w f \\ Zl + \\(E ,B Q )\\ Ll ) 

+ (l+trm\w e+e " 9h Vi +1 d a f \\ + \\Vi +1 d a (E 0l Bo)\\), (3.4) 

where 

a m = ^ + y , t° w > 2a m , if ah > 0, < j < t*? gh . 

Proof. It is divided by the following three steps. For brevity of presentation we 
would only sketch the proof by clarifying how the known techniques in [51 [H] 
can be adopted in the situation considered here. 

Step 1. We claim that there is a time-frequency interactive functional £ lnt (t, k) 
such that 

\£ mt (t,k)\<\f\l 2 + \[E,B}\ 2 , (3.5) 



and 



dt{\f\ 2 



L- 



[E, B}\ 2 + KQ ^£ int (t, k)} + k\{I - P}f\l 

2 Ja^^_| 2 + \b\ 2 + \c\ 2 ) + n\a^a- \ 2 



K \k\ 2 



l+|fc| 2 



(3.7) 



l + \k\ 2 

where kq > is a small constant such that 

|/| 2 2 + \[E,B}\ 2 + Ko ^S' mt (t,k) ~ |/| 2 L2 + \[E,B}\ 2 . 

Compared to the corresponding estimate obtained in [9] , the main improvement in 
(|3.6[) occurs to the coefficient in front of the term \E\ 2 in the dissipation rate. 

Step 2. In this step, we follow the approach in [21] to carry out the velocity 
weighted energy estimates for the pointwise frequency variable. Starting from the 
micro equation 



dt{I - P}f + • k{I - P}f + L{I - P}f 

= -{i-p}(E-OVTm-{i-P} 

one can verity that for £ > 0, 



i£ ■ kPf 



P 



i£-k{I-P}f 



w l {I-P}.f 

< 



L- 



X\k\<\ + K 



1 



l + \k\ 2 



\E\ 



w*{I - P}f 
\k\ 2 



X\k\<l 



1 + lfcl 2 



(a+ + a-,b,c)\ + |a + — a_| 

+ w~ l {I-P}f 



L' 1 



(3.8) 
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In a similar way, starting with the first equation of (|3.I|) . the direct velocity weighted 
energy estimates for the pointwise frequency variable also gives that for t > 0, 



l + \k\ 



w f 



< 



\k\ 



'{I - P}f 



X|fc|>l 



l + lfcl 5 



(a + + a-,b, c)| 2 + |o+ — a_| 2 



l {I-P}f 



L- 



(3.9) 



Now, by choosing K\,K2 > properly small, the linear combination p.6j) + K2 X 
dSHJ) + m x (HU) yields that for £ > 0, 



d t M e (t,k) + nD e (t,k) < 0, 
where Mi(t,k) and Dg(t,k) are given by 

M*(t,fc) = 1 1 /|| i 2 + |[S,B]| 2 + Ko 5R£: int (t,fc) 



(3.10) 



+K2 



w«{I-P}f 



«1 



fc) = ({/-Pi/If, 



^ 1 + ITf 



X|fc|>i> 



(|a + + a _| 2 + |6| 2 + |c| 2 ) + |a + -a_P 



+ 



l + |fc| 2 

1 |£f 



1 + |A;| 2 
Here, notice that for any ^, 



\Ef 



(1 



■|Bp 



«^{/-P}/ > w^{I-P}f 



L 2 



Moreover, set the frequency function p{k) = |fc| 2 /(l + |fc| 2 ) 2 . Then, by consider- 
ing the [1 + ep(fc)i] J -wcightcd estimate on (|3.10[) and also applying the iterative 
technique developed in [6], one has that whenever t > 0, 

M e (t, k) < [1 + ep(fc)t]- J Af £+J+p _ 1 (0, k), (3.11) 

holds true for any t > 0, fc G K 3 , where the parameters p, e and J with 

p > 1, < e < 1, J > 0, CieJ < |. 

are still to be chosen. 
Step 3. For £ > 0, define 

M,(t,fc) = | w £ /l|2 + |[i;,B]| 2 . (3.i2) 

Take a > with m = |a|. We use the splitting 



\k a \ 2 M e (t,k)dk 



k\<l J\k\>l 



\k a \ z M e (t,k) dk. 
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By the estimate p. lip , one has 



|fc|<i 



\k a \ 2 Mt.{t,k) dk 



< 



\k a \ 2 M e {t,k)dk 



< 



k\<l 

2m h i ^„n.\.t\—J 



k <1 



|fc| 2m [l + ep(k)t]- J M e+J+p ^(0, k) dk 



< 



\k\ 2m l + -\k\ 2 t dk sup A/< +J+p _i(0,fc). 
fe|<i L 4 J 



By letting 2J — 2m > 3, i.e. J > m + § = 2cr m , we arrive at 



|fc|<i 



\k a \ 2 Mt(t,k)dk < (1 + (\\w i+J+ ^ 1 hf Zi + 11(^0,50)11 



Here, notice that for any given £ l ° w > 2a m , one can take proper constants J > 2cr„ 
and p > 1 such that £^ ou ' = J + p — 1. Hence, for the low frequency part, 



\k a \ 2 M t {t,k)dk < (i + ty 



|fe|<i 



|(E ,Bo)||ii 



(3.13) 



for ^ 01 " > 2tr m . To estimate the high-frequency part, let j > to be chosen, and 
one has 



\k a \ z M e (t,k) dk 



k\>l 



< 



\k a \ 2 (l + \k\ 2 )M t (t,k) dk 



< 



k\>\ 



\k a \ 2 \k\ 2 M t {t,k)dk 



k\>\ 



< 



\k a \ 2 \k\ 2 [l + ep(k)t}- J M, +J+p ^(0, k) dk 



fe|>i 



which further implies 

\k a \ 2 M e (t,k) dk 
rt 



\k\>i 



< 



< 



k\>l 



1 



4|A| S 



\k a \'\k\ 2 M e+J+p _ 1 (0 1 k)dk 



sup 

|fc|>i 



Notice that 



sup 

|fe|>i 



et 

4\hJ 



\k\ 2j 



et 

+ 4\k\ 2 ~ 



i -J 



\k\ 2 i 



as long as j < J. Thus, we have 
\k a \ 2 Me(t,k) dk 

(i + *)- j (iK+'^-M+^/oll 2 + ||v^ +1 a tt (£;o,i?o)ir) 



fei>i 



< 
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Here, notice that for any given il iah > and < j < one can choose again 

proper constants J > j and p > 1 such that £* 19 ' 1 = J + p — 1. Hence, it follows 
that 



\k a \ 2 M t {t,k) dk 



k >1 



<(! + *)" 



M +1 d a (E Q ,B )\ 



(3.14) 



for £* I9 ' i > and < j < Therefore, in terms of (|3. 12|) . the desired time- 

decay estimate (|3.4j) follows from the combination of (|3.13l) and (|3.14j) . This com- 
pletes the proof of Theorem 13.11 □ 



4. Global a priori estimates 

In this section we will prove Theorem ll.il The key point is to deduce the uniform- 
in-time a priori estimates of solutions to the Vlasov-Maxwell-Landau system 

f fit/ + £ ■ V x / - E ■ ^ 1/2 qi + Lf = S, 

I d t E-V x xB = -(&/ 2 , /+-/_>, r41l 
)d t B + V x xE = 0, 1 ' ' 

{ W x -E= (/i 1/2 , /+-/_>, V»-fl = 0, 

where the nonlinear term is given by 

S = r(/, /) + \q n E ■ U - qo(E + € x B) • V C /. (4.2) 

For that, let (/, E 1 , £?) be a smooth solution to (|4.1[) over the time interval < t < T 
with initial data (/o, Eq, Bo) for < T < oo, and further suppose that (f,E,B) 
satisfies 

< <S 2 , (4.3) 

where is given in (|1.8p and the constant 5 > is sufficiently small. We point out 
that throughout this section, in the case when an undetermined energy functional 
^N,e,\{t) appears on the right-hand side of inequalities, it is always understood to 
take exactly the right-hand expression of (|1.6[) . 

4.1. Macro dissipation. Basing on the previous work [8] and [9], it is a standard 
process to obtain in terms of the following lemma the macro dissipation T>M,mac{t) 
defined by 

V N , mac (t)= J2 l|V^ Q (a±,6,c)|| 2 +|| a+ - a _|| 2 

|a|<iV-l 

The details for the proof of this lemma are omitted for brevity. 

Lemma 4.1. For any integer N with 8 < N < N±, there is an interactive energy 
functional £™ 1 (t) such that 

\£%\t)\< £ (\\d a f\\ 2 + \\d a (E,B)\\ 2 ) 

\a\<N 



GLOBAL SMOOTH DYNAMICS OF A FULLY IONIZED PLASMA 



15 



and 

j t ^\t) + nV N , mac (t) < J2 H^U - P}f\\h + £ N (t)V N (t), 

\a\<N 

forO<t<T. 

4.2. Lyapunov inequality for £^(4). In this section we derive the basic energy 
estimates on £/v/,a(X) m the case of N = N\ and £ = X = 0. Notice that N\ is the 
highest order, and regarding dgf, the weight function takes the form of W—\p\{£) 
independent of time. 

Lemma 4.2. There is an energy functional £isr 1 (t) such that 

— £ Nl (t) + ttD Nl (t) < +t j 1+ # 'DN 1 ,e 1 ,\ (t) +£ Nl (t)£ NoJo - hXo (t), (4.4) 
for < t < T. 

Proof. First of all, it is straightforward to establish the energy identities 

^|(||/|| 2 +\\(E,B)\\ 2 ) + {Lf,f) = (S,f), (4.5) 

and 

Ift E (\\d a f\\ 2 + \\d a (E,B)r)+ (Ld<*f,d"f) 

l<\a\<Ni 1<|«|<JVi 

= J2 (d°S,d a f). (4.6) 

l<|a|<JVi 

By applying the micro projection I — P to the first equation of (|4.1[) . it can be 
rewritten as 

d t {I - P}f + f ■ V x {/ -P}f-E- £m 1/2 + Lf 

= {I-P}S + P£-V x f-£-V x Pf. (4.7) 
After acting dp with |a| + \P\ < Ni and |/3| > 1 to the above equation and fur- 
ther multiplying it by (£) 2 ( 7+2 )l /3 l(9"{/ — P}f, the direct energy estimate gives the 
identity 

±\\8?{I P}f\\fa + {d$Lf, (0 2 ^ +2)l/3| <W - P}f) 

= <^{~e • V X {I -P}f + E- £ M 1/2 }, (0 2(7+2)l/3| ^{/ - P}f) 

+{df 3 {{i - p}s + • v x / - e • v x p/}, (o 2 ^ +2m d^{i - pyf), 

which from Lemma 12.21 implies 
1 H Nl 

2j t ^ Cm E \m i -p}f\\U+ K E \w- p }f\&, w 

m=l |/3|=rn £i|>l 

M+|/9|<ATi |a| + l/3|<JVi 

< j2 \\ d n i - p }f\\i,w\+ E (iiv,3> ± ,Mii 2 + ii3«i?ii 2 ) 

|/3 1 =0 |a|<JVi-l 
\a\ + \0\<Ni 

+ E C ™ E (^{/-p}5,(e) 2 ^ +2 ^i^{j-p}/). (4.8) 

m=l |/3|=m 

|a| + |/3|<JVi 



16 



RENJUN DUAN 



Moreover, from Lemma I4~T1 as well as (|4.3j) . 

±£^(t) + K V Numac (t)< E \\d a {I-P}.f\\l + S 2 V Nl (t). (4.9) 

|a|<JVi 



Then, since S > can be small enough, the proper linear combination of (|4.5[) . 
(|4.6|) . (|4.8j) and (|4.9[) implies that there is an energy functional £N 1 {t) satisfying 
([O]) such that 

(t) + (t) < Z$ (t), (4.10) 

where 

X«(i) = <£,/) + E 

l<|a|<iVi 

+ E C ™ E (^{/-p}5,(o 2(7+2)I ^^{i-p}/). 

m=l |/3|="i 

l"l + l/3|<JVi 

Finally, we claim that 

£ 4f(*)^i(*) + (1+ l 1+ ^ iA,Ao(*) 

+4 / i 2 (*)4 / 1a-i,a (*K 1 2 (*)- (4- 11 ) 

Therefore, f|4.4[) follows from plugging (|4. 1 1[) into (|4. 10[) and applying (|4.3I) and the 
Cauchy-Schwarz inequality. This then completes the proof of Lemma 14. 2 1 □ 



Proof of (|4.11j) . Wc first consider the estimate of X^{t) corresponding to r(/, /) 
in the nonlinear term S. As in [TO], by decomposing / as Pf + {I — P}f and 
using Lemma [2~2l it directly follows that it is bounded up to a generic constant by 

1/2 

£p{ (i)2?/Vi(i)- Recall from the definition of X(t), 

X l / 2 (t) A 
\\V x {E,B)\\ H „ -i< - ■ -lig < 



(l + t)i+* - (i + t)i+* - 
For the zero-order term related to the electromagnetic field, it holds that 

(\qoE ■ U - q (E + £ x B) ■ V f /, /) = (^ £ ■ £/, /) 



< 



\E\.\mPf\ 2 + \{I-P}f\ 2 )dxdt 



<\\E\\-\\(a,b,c)\\ L ~\\(a,b,c)\\ + \\E\\ 



L- 



II \Z\-\{i-P}f\ 2 dxds; 



< e^{t)v Nl {t) + , 1 + t)1+ # v Nu e u x (t). 
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For the d a derivative term related to (E, B) with 1 < |a| < Ni, one has 
(d a (E-Zf),d a f) = C* i (d ai E-£d a - ai f,d a f) 

Oil <a 

< Yl ii^^iu-lieKf) 3 * 2 ^" 01 /!- m^vn 

\a 1 \<N 1 /2.a 1 ^a 

-/+2 



ii^^ii • sup im^d^f \\(o^d a f\\ 

X hp 

oil I >Ni/2 or ai-a 



~ (1 + t)i+« P *iA,Ao(*) + ^(^^A-LAoW^ 2 ^' 

where the Sobolev inequality ||<?||l=° < CH^gH^-i for any function g — g(x) G H 2 
has been used, and we also have used the choice of Ni = |JVo and l\ = \Iq with 
Ao and £q properly large. And in a similar way, it holds that 

(d a {(£xB).v e f},d a f)= J2 c;(Kxa«'j]).v { r"7,37) 

0<Oti <Q 

8 „ / .x _1 /2 / .X ^1 II /.x— 1/2/ 



(1+t)W %A,A n M+^( f )C,-i.J f )^( f )- 



The completely same estimate holds true for the do derivative term related to 
(E,B) with \a\ + \(3\ < Ni and |/3i| > 1, through observing 

{I - P}{\qoE ■ U - q (E + £ x B) ■ V S /} 

= ^<? P • £{I - P}f - 1o(E + (xB). V C {J - P}/ 

+ iq i? • CP/ -»(£ + ?x5). V S P/ 

- P{ig £ • £/ - qo(E + £ x P) • V C /}. 

Therefore, f|4. 1 1[) follows by collecting these estimates. □ 

4.3. Lyapunov inequality for £/Vi.£i.a I n this section we turn to the weighted 
energy estimates on £jVi,£i,A (i)- Notice that due to the regularity doss property of 
the whole system, the weighted high-order energy functional £jVi,£i,Ao(^) can n °t 
be bounded uniformly in time; it has been actually seen from the proof of Lemma 
13.11 for the linearized analysis, cf. (|3.10p . Instead we may expect that £N lt e lt \ (t) 
increases in time. Another trouble arises from the weighted estimate on derivatives 
of the highest order Ni for the linear term E ■ ^/i 1 / 2 . However, these difficulties 
will be resolved by considering the time weighted estimate with the time rate of 
negative power. In addition, although the velocity growth in the nonlinear term 
containing the electromagnetic field could be dealt with through the time-dependent 
exponential factor in the weight function, it is impossible when the electromagnetic 
field gains the differentiation of higher orders since they again could not decay in 
time. Thus, to overcome it, one has to refine the nonlinear estimates in order to 
use the time-decay property of the lower-order energy functional £N ,e -i.\ a with 
the higher-order algebraic velocity weight. 
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Lemma 4.3. There is an energy functional £isf 1 ,£ 1 ,\ (t) such that 

^£jViA,Ao(*) + K ' D N 1 ,l 1 ,\ {t) < E Nl ,i 1 ,\ Q {t)E Noi i -i^\ Q {t) 

+ (d a E.^ 2 ,w\, Xo d a f), (4.12) 
H=iVi 

for < t < T. 

Proof. Starting from the first equation of (|4.ip . the energy estimate on d a f with 1 < 
M < Ni weighted by the time-velocity dependent function w_i lXo = ^ (f, £) 
gives 

~ E + E (urf,^*,**/) 

l<\a\<N ± l<l«|<JVi 

+ 7 I ^T? IKWH 2 -. 1 ,Ao= E (d a S,w 2 _ iuXa d»f) 

l<|o|<JVi 

+ ^ {d a E-^'\wl tlM d a f). (4.13) 

1<|q|<JVi 

Similarly, from (|4.7p . one has the weighted energy estimate on {/ — P}f 
P}f\W,x + (Lf,w 2 _ iuXo {I P}f) 



9X0 ^\\{0{I-P}f\\-i u x 



= <{/ - P}S, w\ M {I - P}f) + (E ■ &V 2 ,w 2 _ iuXa {I - P}f) 

+ (P£ • V„/ - f ■ V a P/, u£* llAo {J - P}f), (4.14) 

and the weighted energy estimate on {I — P}dpf with \a\ + \j3\ < N\ and \/3\ > 1 
1 d Nl 

2j t ^ Cm £ w-n/iijw,*, 

m=l |/9|=m 

|a|+|0|<JVi 

+« E (ll W - P }/lln^|-, 1> A + (1+ A ° 1+ J I(Ogg{/ - P}/||) 

M+|;8j<A7i 

< J2 W d n i - p }f\\nM-^Xo+ E (l|V^(a ± ,&,c)|| 2 + ||^|| 2 ) 

101=0 |a|<JVi-l 
|a|+|0|<JVi 

AT! 

+ E C " E W-P}S,^- tu}i0 ^{I-P}f). (4-15) 

m=l |/3|=m 

|a| + |0|<JVi 
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Then, the proper linear combination of (|4.5jl . (|4~6]) . (|4.13j) . (|4.14[) and ()4.15|) implies 
that there is an energy functional £ n 1 ,i 1 ,\ (t) satisfying ()1.6|) such that 

+ {d a E-^ l '\wl tlM d a f). (4.16) 

M=JVi 

where 

l<|a|<]Vi 

+ ({/-P}5,^ 1)Ao {J-P}/)+ E (^wl^ff"/) 

l<|a|<ATi 

Ni 

+ E C » E (^u-m^i-^^w-n/)- (4.i7) 

m=l |/3|=to 

|a| + |/9|<JVi 

We now claim that 

T N ) 1 ,e 1 ^ ( t ) < ^oA.AqW^Wi.AoC*) +^A,Ao^) f M ) 2 £o-l,Ao^) 23 W 1 2 £ 1 ,A ( < ) 

+ ^(l + 1+ "l|V :c ( J B,i?)|| ff «o- 1 ^ 1 A,A (i)- (4-18) 

Ao 

Notice that the first and third terms on the right are bounded up to a generic 
constant by (WiViA.AoW due to the definition of X(i) and the assumption X(t) < 
S 2 . This hence simplifies (|4.18[) as 

4? A ,a (*) Z ^Ka.Ac W4!a-i,a W^a.Ao (*) + ^A,A (i)- (4-19) 

Therefore, since ^ > is small enough, by applying the Cauchy-Schwarz inequality 
to the first term on the right-hand side of (|4.19p and then putting it into (|4.16l) . 
(|4.12[) follows. This completes the proof of Lemma [4.31 □ 

Proof of (|4.18j) . For brevity, we only present the estimate of the fourth term on 
the right-hand side of (|4.17|) since the estimate on other terms is simpler or follows 
in the completely same way. Take a with 1 < |a| < Ni. For the inner product 
term related to d a T(f,f), by using (|2.2[) in Lemma [2.2[ considering the cases of 
\a'\ < Ni/2 and \a'\ > iVi/2, applying the Sobolev inequality \\g\\ L °° < C\\\7 x g\\ H i 
to the lower-order derivatives, and noticing the choice of Ni = §iVo and £\ = ^£q 
with A^o and £q properly large, it follows that 

(d a r(fj),w 2 _t uXo d a f) 

~ {^A ,«o,A (*) I? Afi,fi,Ao( t ) +^iVoA-l,A 
~^~{^N 1 ,e 1 ,X (^)^Afo,fo-l,A 
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Next, for the term E ■ £/ in S, one has 



< J2 W d °" E h° 



|ai|<ATi/2 



|^|^ 1)Ao (|^-«vr + l^/| 2 )^ 



|ai|>ATi/2 
or Qi-a 



1 



<^(l + i) 1+ *||V a (£?, BJII^o-i^ AiAo (t) 
Ao 

(«a,a W- 

For the term (E + £ x _B) ■ Vf / in 5, the difference point is that it contains the 
velocity derivative of order one. One can deduc that 

(d a [(E + txB)-V c f],w 2 _ eiiXo d a f) 
= ((E + txB)-V s w 2 _ eiM ,-±\d a f\ 2 ) 

+ E C^((d^E + ^xd^B)-V 6 d a -^f,w 2 _ iuXo d a f), (4.20) 

0<ai <a 

Here, it is straightforward to see that the first term on the right is bounded in a 
rough way by 

(0 



< J_(l +t )l+*llv a£ ;ll BI 

Ao 



Ao(Q 2 
(1 + 



< ( 1 + t) i+*\\V xE \\ H1 V NlA , Xo (t). 
Ao 



The second term on the right-hand side of (|4.20j) can be estimated as follows. When 
| ofi j < N±/2 and ot\ ^= a, it is bounded by 

C\\ff"{E,B)\\ L - jj (0 1 - h+2 H^-, 1 ^ed a - a \f\ 2 +w\ M \d a f\ 2 )dxd^ 

J Jl 3 x E 3 

< \\V x d a -(E,B)\\ H i f f (O^^i^.AoV^"-" 1 /! 2 + \w- tl , Xo d a f\ 2 )dxd£ 

< ^-(l + tJ^IIVx^.flJH^o-iXJjVL/LAoW. 
Ao 

where we have used —3 < 7 < —2, and when > N\/2 or a% = a, it is bounded 

by 

C\\d ai (E,B)\\ -sup (O^^^.AoVe^-" 1 / I|<0^-^a <9 q /II 

a; ^ 

Collecting all the above estimates, (|4. 18[) holds true for the fourth inner product 
term on the right-hand side of ()4.17j) . This proves (|4.18|) . □ 
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Now, we are ready to obtain the closed estimate on the first portion of the 
time- weighted energy norm X(t) in the following 

Lemma 4.4. It holds that 

sup {£ Nl (s) + (1 + sy^£ NueiM {s)\ + f V Nl (s)ds<Y 2 + X 2 {t), (4.21) 

0<s<t L ' Jo 

forO<t<T. 

Proof. In fact, the time integration of (|4.4j) gives 

£ Nl (t)+ f V Nl (s)ds<Y 2 + S f (l + 8)- 1 -*V Nl j u x (s)d8 
Jo Jo 

+ [ £ Nl (s)£ NoA) _ 1M (s)ds. (4.22) 
Jo 

Furthermore, from multiplying f|4.4[) by (1 + t)~ t0 and then taking the time inte- 
gration, it follows that 

(l + t)- £0 £ Nl {t)+ f (1 + s)- f -°V Nl (s) ds + f (1 + s)- 1 - eo £ Nl ( S )ds 
Jo Jo 

<Y 2 + 6 f (l + sy^^V^^^ds 
Jo 

+ [ (l + s)- e °£ Nl (s)£ No ,i - 1 , Xo (s)ds. (4.23) 
Jo 

Combining (|4~2"2"j) and (g^) gives 

£ Nl (t)+ f V Nl { S )ds+ f (l + S )- 1 - eo £ A r 1 (s)ds 

Jo Jo 
<Y 2 +sf {l + s)~ 1 -< > V NueiM (s)d S + f £ Nl {s)£ No ^ hXo (s)ds 



10 JO 

' Nl ,i lt \ (s) ds, (4.24) 
where to obtain the second inequality, we have used 



<Y 2 +X 2 (t) + S f (1 + s)" 1 -^ 
Jo 



sup {£ Nl (s) + (1 + s)z£ No . eo - 1M (s)} < X(t). 

0<s<t 

From (|4.12p . multiplying it by (l + t)~( 1+e °)/ 2 and taking the time integration yields 

(1 + t) —£ Nl .i 1 .x {t)+ / (l + s) —V NlMM (s)ds 

Jo 

+ / (1 + sy^£ Nl , eiAo {s) ds < F 2 + / (1 + s)-^ £ Nl4liKo ( s )£ No}io _ liXo (s) ds 



Jo 

t 



+ E / (l + s)-^(d a E-^ 2 ,w 2 _ luXo d a f). (4.25) 

The second term on the right is bounded by X 2 (t) by noticing again from the 
definition of X(t) 

l + ep 3 

sup {(1 + s) 2 £n u h u \ {s) + (1 + s)2£ No ,e -i,\ (s)} < X(t). 

0<s<t 
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By the Cauchy-Schwarz inequality, the right-hand third term of (|4.25[) is bounded 
up to a generic constant by 

E [ (l + s)- 1 - eo \\d a E\\ 2 + \\(0 1 ^d a f\\ 2 ds 

< f (l + s)- 1 - eo £ Nl (s)ds+ f V Nl (s)ds. 
Jo Jo 

Then, in terms of the above estimates, taking the sum of (|4.24[) and f|4. 25[) and 

using the fact that 5 > is small enough, we arrive at 

£*ri(t) + (1 + t)~^ £jviA,Ao(*) + / VNi(s)ds 

Jo 

rt rt 

+ / {l + s)-^V NlMM {s)ds+ / (l+s)- 1 - £0 ^ 1 (s)ds 
Jo Jo 

+ / (l + s)-^£ NuiuXo (s)ds<Y£+X 2 (t). (4.26) 
Jo 

Therefore, (|4.21[) follows, and then this completes the proof of Lemma T4. 41 □ 

4.4. Lyapunov inequality for £N a ,e ,\ a (t)- In this section we devote ourselves to 
obtaining the Lyapunov inequality for £/v ,£ 0i a (0- It is simpler compared to deal 
with the high-order energy functional £jvi,£i,Ao(*)i because the derivatives of the 
electromagnetic field of order up to iVo decays in time by the definition of X(t). 

Lemma 4.5. For any £ with £q~1 < £ < £q, there is an energy functional £n ,i,\ (t) 
such that 

|a|=JV 

for < t < T. 

Proof. From the completely same procedure to obtain the energy inequality (|4.16[) 
for £/vi,£i,A {t), one can also verify that for I with £q — 1 < £ < £g, there is an energy 
functional £N a .i,\ (t) satisfying (|1.6[) such that 

^oAA (*) + «»iVoAA„(t)<^ 1 Ao(*)+ E (9 a E.^ 2 ,w 2 _ lXo d a f). (4.28) 

|a|=JV 

where 

<W*) = <£,/>+ E (d a S,d a f) 

l<\a\<N 

+ ({I-P}S,w\ Xa {I-P}f)+ E (d a S,w 2 _ eM d a f) 

l<\a\<N 

N 

+ E C - E {d${I-P}S,w 2 m _, M d a p {I-P}f). (4.29) 

m=l \j3\=m 

\a\ + \/3\<N 

(2) ______ 

Compared to the estimate on the similar functional X N $ x (t) given by (|4.17p in 
Lemma 14.31 the estimate on X No t Xq (t) above becomes easier due to the bounded- 
ness of £n ,i,\ {^) uniformly in time and also the time-decay of derivatives of the 
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electromagnetic field up to order Nq. In fact, one can claim that 

+ y o (l + t) 1+ ' > \\V x (E,B)\\ HNo - 1 V NoAXo (t). (4.30) 

Noticing from the definition of X(t) 

sup {£ NMo (s) + (1 + S ) 2 ( 1+ ^||V :c (^ )J B)||^ _ 1 } < X(t) < <5 2 , 

0<s<t 

it further follows that t Ao (t) is bounded up to a generic constant by <52?jv .£.A (*) ■ 
By putting this estimate into (|4.28p and applying the Cauchy-Schwarz inequality 
to the right-hand second term of (|4.28[) as 

H=iVo \a\=N " \a\=N 

for -q > small enough, it follows that 

T £N ,tAo(t) + KV NoAXo (t) < (S + v)T> NoA x (t) + - W daE W 2 > 
at r) , , — ' 

which thus implies (|4.27p since 6 > and rj > can be small enough. This completes 
the proof of Lemma 14.51 □ 

________ ______ 

Proof of (|4.30p . It is similar to the proof of (|4.18[) for A (t). For those inner 

product terms from T(f,f) in (|4.29[) . one can directly apply Lemma \2. 2 1 to verify 

1 /2 

that they are bounded by C£^ ( A (t)2? Ar ,^,A W ■ F° r those terms related to (E, B) 
in (|4.29p . from the completely same process as for dealing in the proof of (|4.18p 
with the case when (E, B) gains the differentiation of lower-order, it follows that 
they are bounded by C^(l + t) 1+§ \\ V X (E, B)\\ h n -iV NoAXo (*)• Therefore, (jOO)) 
is proved. □ 

4.5. Decay of electromagnetic fields and macro components. In this step, 
we will use directly the Duhamel's principle to obtain the time-decay of the electro- 
magnetic field (E, B) and the macro components (a, b, c) up to the low-order Nq in 
terms of the time-decay of the weighted high-order energy function £N 1 -3,e 1 -i,\ a (t) 
which follows from the boundedness of X(t). 

Lemma 4.6. It holds that 

sup {(l + S )^\V x (E,B)\\% No ^ + (l + s)i\\(a,b,c,E,B)f} <Y 2 + X 2 (t), (4.31) 

0<s<t 

for < t < T. 

Proof. Recall the mild form 

U(t) = A(t)U + [ A(i- s)[S(s), 0,0] da, (4.32) 
Jo 

which denotes the solutions to the Cauchy problem on the Vlasov-Maxwcll-Landau 
system (|4.ip with initial data t/o = (/o, Eq, Bq), where the nonlinear term S is 
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given by (|4.2j) . The linearized analysis for the homogeneous system in Theorem l3.ll 
implies 

\\VxPE,B{Mt)Uo}\\ H "o-i < (l+t)-i(\\w^ w fo\\ Zl +\\(E ,B )\\ L i) 

+ (l + f)~* (lk 4<9 "v|a a / || + \\\/ld a (E a ,B )\\), 

l<\a\<N 

where Pe.b means the projection along the electro and magnetic components in 
the solution (f,E,B), w = w(£) is dehned by (|3.3p . and constants £ l ^ w , l^ 9h are 
chosen to satisfy 

and also t l £ w , £^ l9h are sufficiently close to 3/2 + Nq and 5/4, respectively. By 
interpolation of derivatives, 

||V a Pte,B{A(t)^o}|| ff *o-i < (l+t)-i(\\w^ W fo\\ Zl + \\(Eo,Bo)\\ L i) 

+ + (\\w e " 9h d a M\ + \\d a (E Q ,B )\\). 

3<|a|<JV +3 

Applying this time-decay property to the mild form (j4.32[) gives 

rt 

\\V x (E,B)\\ H * -i<(l + t)-lY + / (l+t- s )-i\\w el ™S( S )\\ Zl ds 

Jo 

+ / (l + t-s)-i \\w ih 3 t9h d a S(s)\\ds, (4.33) 

3<|a|<JV +3 

where we have used w(£) = w_i(£) and the definition (|1.9[) for Y . As in [7], it is 
straightforward to obtain 

\\v&"S{t)Ui + E \\^ i9h d a S(t)\\ <£ Nl . 3A . 1M (t). 

3<\a\<N +3 

Here, we have used the choice of N±, l\ by N\ = §Ao, i\ — \io with No and to 
properly large. Recall X(t) norm, and hence 

fjvi-sA-i.AoOO < (1 + s)-^X(t), 0<s<t. 
Plugging these estimates into (|4.33[) . the further computations yield 

sup {(l + S ^\\V x (E,B)\\ 2 HNo ^}<Y 2 +X 2 (t). (4.34) 

0<s<t 

Moreover, to obtain the time-decay of || (a, 6, c, E 7 B)\\ , we use the linearized time- 
decay property 

\\P f {A(t)U }\\ + \\P E>B {A(t)U }\\ < (l + i)-l(||^r /o || Zi + \\(Eo,B )\\ L i) 

+ (1 +t)-i(\\w i * i9 \lf \\ + \\wl(Eo,B )\\), 

where Pf means the projection along the /-component in the solution (f,E,B), 
and constants t l £ w , if 9 * 1 are chosen to satisfy i l £ w > 3/2, l h ^ h > 3/4 and also £ l £ w , 
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^ gh are sufficiently close to 3/2 and 3/4, respectively. Therefore, in the completely 
same way for estimating || \7 X (E, B)\\ h n -i in (|4.34j) . one has 

sup {(1 + s)?\\{a 1 b,c,E,B)\\ 2 }<Y 2 +X 2 (t). (4.35) 

0<s<t 

Thus, combining (|4.34p and (|4.35[) gives the desired estimate (|4.3f|) . This then 
completes the proof of Lemma 14.61 □ 



4.6. Bound of £N .io-X (t) an£ i decay of £tv ,£o-i,AoW- Basing on those estimates 
in the previous two sections, we can use the time-weighted estimate together with 
an iterative trick to obtain in the following lemma the boudedness of £N .e a .x (t) 
and also the time-decay of £N ,e -i,\ (t). Notice that loss of one in the velocity 
weight index of £at ,£ -i,Ao (t) results essentially from the long-range degenerate 
property of soft potentials for the Landau operator. 

Lemma 4.7. It holds that 

sup {£ NMo (s) + a + s)^N ,e -i.x (s)} <Y a 2 + X 2 (t), (4.36) 

0<s<t 

for < t < T. 

Proof. Recall the Lyapunov inequality (|4.27|) for £N ,t,x (t) with £q — 1 < £ < £q. 
First of all, from the time integration of (|4.27p with I = £q, 

£*Mo,Ao(*)+ [ t V NaAlM (s)ds<Y 2 + ]T f\\d a E\\ 2 d S . (4.37) 

J ° \ a \=N Jo 

Due to Lemma l4.6[ the second term on the right is bounded up to a generic constant 
by Yq + X 2 (t), and so is £ n ,£ ,^o ( s ) for < s < i. The rest is to estimate the 
time- weighted part on the left-hand side of (|4.36|) . For that, let e > be fixed small 
enough. From multiplying (|4.27p with £ = £q — 1/2 by (1 + t)^ +e and then taking 
the time integration, 



(l + t)^£ NQAl _ hM {t)+ / {l + S )^V NoAs _ hM {s)ds<Y 2 

J 

+ [\l + s)-^£ Notto _ h , Xo (s)ds+ [\l + s) h+e \\d a Efds, (4.38) 

and in a similar way, starting from (|4.27|) with £ = £q — 1, 



(l + t)2+ £ £ JWo _ 1 , Ao (t)+ / (l + s )^V No ^ 1M ( s )ds<Y 2 

Jo 

+ [ (l + s)i +e £ No , eo - 1M ( S )d s + V [ (l + s)^\\d a E\\ 2 d s . (4.39) 

JO ,_, »r JO 



\a\=N ' 
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Using the relation between the energy functional £N,e,\(t) and its dissipation rate 
1^N,l,\(t), the proper linear combination of (|4.37|) . (|4.38[) and (|4.39|) yields 



(1 + t)* +e fjv oA -i,Ao(*) + / V NMa {s) + (1 + s)* +t V Na ^-i M {s)ds 

Jo 

<Y 2 + f\l + S )^\\(a,b,c,B)\\ 2 d S 
Jo 

+ V / (l + s)i +i \\d a (E,B)\\ 2 ds. (4.40) 



\a\=N 

Again from Lemma 14.01 one has 



[\l + s)? +e \\(a,b,c,B)\\ 2 ds + f t ^ + ^ +f -\\d a (E,B)\\ 2 ds 

J ° \ a \=N J ° 

< f\i + ds [Y 2 + X 2 (t)} < (1 + ty\Y 2 + X 2 (t)}. 
Jo 

Using this, it follows from (|4.40j) that 

sup {(1 + s)^ No j . 1Ao (s)} < Y 2 + X 2 {t). 

0<s<t 

Therefore (|4.36|) holds true. This completes the proof of Lemma 14.71 □ 

4.7. Bound of £ N±-i,i±,\ m this section we obtain the uniform-in-time bound- 
edness of the energy functional £n 1 -i,£ 1 ,\ (t)- Notice that this is consistent with 
(13.101) in the linearized analysis. The main observation in the nonlinear analysis is 
that those remaining terms in the energy inequalities are time-space integrable. 

Lemma 4.8. It holds that 

sup Sn 1 -iaM 8 )+ f V Nl - h e u x {s)ds<Y 2 +X 2 (t), (4.41) 

0<s<t Jo 

for < t < T. 

Proof. Similarly for obtaining (|4.16|) . one has 

—£ Nl -lj lM (t) +/cDjv 1 -1,4,A (*) ^^Vi-l^AoW 

+ J2 (d a E-^ 2 ,w 2 _ iuXo d<*f), (4.42) 

|a|=JVi-l 

where I Nl _ x tl Xo (t) is defined by (|4.17[) with N\ replaced by Ni — 1. The first term 
on the right can be bounded by 

+ ^(1 +t) 1+ »\\V x (E,B)\\ HNo -iV Nl - 14l , Xo (t) 

*0 

+ Ci(<Vu„W<ii 1 ,A„W- (4-43) 
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Further using 

sup {£ Nl -i(s) +S Nl -i AuXo {a) + (1 + s) 2 ^\\V :c (E,B)\\ 2 HN ^ 1 

0<s<t 

+ (l + s^£ NoA ,. liXo (s)}<X(t)<S 2 , 

it follows that 

^S-iA.AoW £ 5V Nl ^ eiM (t) + (l+t)-iX(t)V^_ 1>iuXo (t). (4.44) 

From the Cauchy-Schwarz inequality, the right-hand second term of (|4.42[) is esti- 
mated by 

|a|=JVi — 1 | a | = JVi — 1 ? ' 

< r)-D Nl - 1AM (t) + -V Nl (t), (4.45) 
V 

for any r\ > 0. Then, by applying again the Cauchy-Schwarz inequality with 77 to 
the right-hand second term of (|4.44[) , plugging the resultant estimate together with 
(|4.45p into (|4.42[) , and choosing 77 > small enough, one has 

j t £ Nl . 1MM {t) + fs2?Ar 1 -iA,A (<) £ fjvi(*) + (1 + t)-ix 2 (t). (4.46) 
Recall that from P~2"6"|) . 

ft 

V Nl {s)ds<Y 2 +X 2 {t). 



Therefore, (|4.41[) follows by the time integration of (|4.46[) . This completes the proof 
of Lemma 14.81 □ 



4.8. Decay of £N 1 -3,£ 1 -i,\ a {t) an d £^-2^). To obtain the closed estimate on 
the energy norm X(t), it remains to obtain the time-decay of the high-order energy 
functional £N 1 -3.i 1 -i,\ (t) and £pj 1 _2(t) through the time-weighted estimate as 
well as the iterative trick as for dealing with £jv 0j £ _i,A (i) in Lemma l4~7l Notice 
that loss of three and loss of one in the smoothness and velocity weight indices in 
£n 1 -3,£ 1 -i,x (t) respectively result from the regularity-loss of the electromagnetic 
held and the degeneration of collisional kernels for soft potentials. 

Lemma 4.9. It holds that 

sup {(1 + 8)i[£ Nl - 3i e 1 -i,M + Z Y o + * 2 W> ( 4 - 47 ) 

0<s<t 

for < t < T. 

Proof. First recall from Lemma 14.41 and Lemma |4. 



£N 1 (t)+£ Nl -i.e 1 .\ {t)+ f V Nl {s) +V Nl - 1 .e 1 . Xo (s)ds <Y 2 + X 2 (t). (4.48) 

Jo 

To obtain the time-decay of £jVi-3,.«i-i,;\o W and £j^ 1 _ 2 (t), we will make the time- 
weighted estimate. For brevity of presentation we write 

4?L(*)= E {d a E-^' 2 ,w\ Xo d a f). 

\a\=N 
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From the proof of Lemma 14.21 and Lemma 14. 3[ cf. ()4.10|) and (|4.16[) , one has the 
Lyapunov inequalities 

f Jt £Nl - l{t) + kVn ^ 

Those terms on the right can be estimated as follows. Similar to (|4.11j) . it holds 
that 

+ ^ 1 -l(*)45 -l ) A (*)^ 1 -lW- 

Here, noticing that £^/ i 2 _ 1 (i) < X l l 2 (t) < S is small enough for the first term on the 
right and applying the Cauchy-Schwarz inequality to the third term on the right, 
it then follows from the first equation of (|4.49[) that 

at 

(i). (4.50) 

Moreover, similar to (|4.43[) . it holds that 

+ f (l + 1+,5 ||V ;r (i?,i?)|| ff »o- 1 I? J v 1 - 2A -i,A W 

+ ^JVi-2(*)^JV ,^ -l,Ao( i ) I) JV 1 -2A-i,Ao^' 

which by using -X'(t) < <5 2 for the first two terms on the right and the Cauchy- 
Schwarz inequality for the last term, further implies 



1 



(2) 



1 



~ (H# ffl _ 2A i, Ao (t) + -fM-a(*)^oA-i.Ao(*)> (4-51) 
for 77 > 0. Again from the Cauchy-Schwarz inequality with 77 > 0, 

fSl-W-i*,® $ E (vUO^dVf + kd^Ef). (4.52) 

|a|=JVi-2 

Then, by plugging (|4.51[) and (|4.52[) into the second equation of (|4.49[) . taking the 
sum of the resultant inequality multiplied by a proper small constant K3 > and 
another inequality (|4.50[) . and using smallness of 5 > and 77 > 0, one has 

+ WS^-Vi-I.AoW} + JVi-l(t) + «3^iV 1 -2A-i,Ao( i )} 

~ (1 + m+^ - PjVl ~ 1 ^ 1 ' A °^ +^JVi-i(*)^JV ,^ -i,Ao(*)- (4.53) 
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Further from multiplying it by (1 + t)^ +e with e > fixed small enough and taking 
the time integration, it follows 

(l + tjH^^t) + £ Ni _ 2M _ hXo (t)} 

+ J Q ( 1 + s f 2+e {V Nl -i(s) + V Ni _ 2/i _ hM {s)} ds 

<Y 2 + f 5(1 + s)-i-^V Nl ^, euXo ( s )ds 
Jo 

+ / (1 + s)^ +c £ Nl ^i(s)£ Not £ - 1 .x (s)ds 
Jo 

+ [ (l + S )-i +e {S Nl -i( S )+S Nl -2M-i^( s )} ds - ( 4 - 54 ) 



Here, since e > is small enough, the second term on the right is bounded by 
Yq + X 2 (t) directly by (|4.48[) , the third term on the right is bounded by 



S 2 sup {(l + s)^£ Nl ^( S )}, 

0<s<t 



due to the fact that 



sup {(l + s)*£ No ,i -i,x (s)} < X(t) < 5 2 , 

0<s<t 

and the fourth term on the right is bounded by Y 2 + X 2 (t) by noticing 

£m-l(t) +£ Nl _ 2/l _ h Xo (t) <V Nl (t) +£ JVl _ 1/l , Ao (i) + \\(a,b,c,B)\\ 2 , 
and further using (|4.48|) as well as Lemma l4~6l Hence, we arrive from (|4.54l) at 
sup {(l + s)^^.^^^ , ( S )]} 

0<s<t 2 

+ I (1 + a)* +e {ZWi-i(«) + $,a»} ds Z Y o 2 + X 2 (t). (4.55) 

In a similar way to obtain (|4.53[) . starting with the Lyapunov inequalities 



d 

'dt" 

-( 2 ) M -L 



-j;£N 1 -3,e 1 -l,\ (t) + «I?Ar 1 _3^ 1 _i i A (t) 



one can prove 

^l{^iVi-2(*) + ^^Wi-S^i-l.Ao^)} + n{T>N 1 -2{t) + K i VN 1 -2,e 1 -l.\ (t)} 

— ^JV!-2A-i,Ao(*) + £jVi-2(t)£/VoA>-l,Ao(*)i 



(1 + *) 
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for a properly chosen constant K4 > 0. Further multiplying it by (1 + t) 2+e and 
taking the time integration gives 

(l + ^i+^JV^W + ^-sA-i.AoC*)} 

+ / (I + s)§ +e {V Nl - 2 (s) +V Nl -3, tl -i,\ (s)}ds 



Jo 

< Yo + J* 80- + s) i+e -*Vm-2M-hM( s ) ds 
t 

(1 + s)2 +c £ Nl _ 2 (s)£ Nai i - 1 ,\ a (s) ds 







+ / (1 + s)* +e {£ Nl _ 2 (s) + e Nl - attl -i M (s)} ds. (4.56) 



Here, notice again that e > is a fixed constant small enough. Then, the second 
term on the right is bounded by Y~ 2 + X 2 (t) by (|4.55p . the third term on the right 
is bounded by X 2 (t) due to 

sup {(1 + s)i[£ Nl - 2 {s) + £N ,e -h\o( s )}} <X(t), 

0<s<t 

and as before, the fourth term on the right is bounded by 

c(i + ty[Y 2 + x 2 (t)} 

by noticing 

£ Nl -2(t) +£ Nl - 3 , il - 1 ,x (t) < V Nl .x{t) +D JVl _ Vl _i iAo (t) + ||(a,6,c,B)|| 2 , 

and further using (|4.55|) as well as Lemma 14.61 Therefore, the desired inequality 
(|4.47|) follows by putting these estimates into (|4.56[) . This then completes the proof 



of Lemma [491 □ 

4.9. Global existence. We are now in a position to complete the 

Proof of Theorem ] Recall X-norm. From Lemma f4"4l Lemma Lemma |4~71 
Lemma HH] and Lemma l4~9l it follows that 



X(t) < Y 2 +X 2 {t). 

Since Yo is sufficiently small, (| 1 . 10[) holds true. The global existence follows. □ 
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